Evolving networks with disadvantaged long-range connections 
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Abstract. We consider a growing network, whose growth algorithm is based on the preferential 
attachment typical for scale-free constructions, but where the long-range bonds are disadvantaged. 
Thus, the probability to get connected to a site at distance d is proportional to d~ a , where a is a 
tunable parameter of the model. We show that the properties of the networks grown with a < 1 are 
close to those of the genuine scale-free construction, while for a > 1 the structure of the network is 
vastly different. Thus, in this regime, the node degree distribution is no more a power law, and it 
is well-represented by a stretched exponential. On the other hand, the small-world property of the 
growing networks is preserved at all values of a. 

PACS numbers: 89.75.-k, 05.50.+q, 89.75.Hc 



Complex weblike structures (the small-world or scale- 
free networks) have recently become an object of exten- 
sive investigation, and in the last years a great success 
in understanding the properties of these structures was 
achieved (see Ref. ||] as a review). Apart from appeal- 
ing mathematics, this recent interest is due to the fact 
that many natural and technological systems, like poh 
mer networks , the science collaboration network | 
or networks of chemical reactions in a living cell 
seem to be organized according to some internal princi- 
ples. Thus, the Internet ||, the network of human sexual 
contacts or the WWW [jll) possess a similar struc- 
ture, e.g. are they all based on the preferential attach- 
ment of the newly introduced nodes to the highly con- 
nected old ones. All these networks show the small- world 
property: the typical distance (in terms of the number 
of intermediate connections) between two nodes grows 
logarithmically with the web's size. 

One of the prominent examples of a mathematical 
model of such a growing network is the scale-free (SF) 
construction of Barabasi and Albert fj],[l2j ; and one of its 
most interesting properties is the very specific form of the 
probability distribution of the degree of nodes (i.e. of the 
numbers of bonds connecting any given node i with other 
ones in the network): P(k) oc k~ 3 [jjjl^-jla] . Many mod- 
els have been presented, based on the same two most 
important ingredients: growth and preferential attach- 
ment. Examples are models with an accelerated growth 
of the network Jl6,|i7| , models with a nonlinear prefer- 
ential attachment |15| , with nodes provided by a initial 
attractiveness |13|,|18|, with growth constraints as aging 



and cost |19 20(1, models that have a competitive aspect 
of the nodes |2l| ] , or models of networks that incorporate 
local events as the addition, rewiring or removal of nodes 
or edges p2| . 

The SF-construction may be a reasonable approxima- 
tion for such world-spanning networks like one of the In- 
ternet's information transmission channels or one of the 
formal links of WWW. On the other hand, in many sit- 



uations (like in a network of human sexual contacts) a 
connection means a physical contact, i.e. means that 
the contacting individuals, representing the nodes of the 
network, have to occur at the same site and at the same 
time, thus introducing a clear geographical aspect. In 
what follows we present a simple model taking into ac- 
count this metrical ("geographical") aspect, where the 
probability to connect two nodes depends both on the 
number of connections that the nodes already have (as 
in the genuine SF-construction), and on the distance be- 
tween them. That is, we treat an emerging network in a 
metric space. In this emerging network the probability 
that a newly introduced node n is connected to a pre- 
viously existing node i is proportional to the number fcj 
of the already existing connections of node i (preferen- 
tial attachment prescription), but on the other hand the 
too long bonds are disadvantaged, because this probabil- 
ity depends on the Euclidean distance di n between the 
nodes n and i as d~f*, (clearly, a "scale- free" function), 
with a > 0. 

Based on extensive numerical simulations of a one- 
dimensional situation, we show that even if the length 
penalties are mild, the model exhibits properties which 
differ strongly from those of the usual scale-free net- 
works. Thus, the corresponding degree distribution func- 
tion P(k) depends strongly on a. We show, in par- 
ticular, that for a < 1 the behavior of P(k) is simi- 
lar to the behavior of the SF model without penalties, 
so that asymptotically P{k) oc k~ 3 , (which distribution 
possesses a mean, but no dispersion, and corresponds 
to strong, universal fluctuations). On the other hand, 
for a > 1 the behavior of P(k) is well-described by a 
stretched-exponential P(k) oc exp(— bk 1 ), with the power 
7 depending on a, so that the fluctuations in k are rather 
weak. We discuss the reasons for such a dramatic change, 
being rooted in the probability of connection between the 
nodes as function of the distance, and the overall struc- 
ture of the emerging network, preserving its small-world 
nature even at large (probably at all) a- values. 
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FIG. 1. Networks generated using the simulation prescription, Eq.(2), with different values of a: (a) a = 0, (b) a — 1.5 
and (c) a = 15. All three examples have 300 edges, L = 10 6 , N = 102, and m = 3. Note the change in the appearence of 
the networks. The network (a) is a genuine SF construction while (c) strongly resembles the Watts and Strogatz's small-world 
network. 



We start from a one-dimensional lattice of L sites, 
spaced by a unit distance and apply cyclic boundary con- 
ditions. On this structure we will let our network grow, 
so that each lattice site will be a possible location of a 
network's node. We denote by t%i the position in the lat- 
tice of a node i. The distance dij between any two nodes 
i and j is defined as: 

d i3 = min{\ni - rij\, (L - |n* - rij\)}. (1) 

Let us now construct the network. First, we choose 
randomly an even number too of sites from the lattice 
and we bind them in pairs with one bond each. This will 
be our initial condition. That is, at t = 0, our network 
will consist from Too nodes connected in pairs. As in the 
SF model we will add at every time step a new node to 
our network (linear growth). We proceed according to 
the following rule: at every time step we choose at ran- 
dom a free site of our lattice, and pose the new node 
there. This new node is then connected through m edges 
(to < mo) with to different nodes already present in the 
network. After t time steps the algorithm results in a 
network with i + TOo nodes and mt + mo/2 edges. In con- 
trast with the SF model, the probability II for the new 
node n to be connected to an old one i will depends not 
only on the number of edges fcj which i already possesses, 
but also on the distance di n between them: 

n(k i ,d in ,a)= k \ d ™ (2) 

l^j K 3 ' a jn 

Here the sum in the denominator goes over all nodes 
in the system except the newly introduced one and a 
is a real non-negative parameter describing the distance 



penalties. For large a, the probability of connection be- 
tween two distant nodes is very small. On the other hand, 
for a very small a the probability is almost independent 
from the distance. In the case a — our model reduces 
to the genuine scale-free one. Note that our model is to 
some extent also scale-free: the connection probabilities 
depend only on the relative distances. 

Our initial condition is slightly different from one of 
Barabasi and Albert, where the initial too nodes are not 
connected: in our case all nodes introduced at t — have 
exactly one edge, which allows to use Eq.(|2|) from the 
very beginning. This simplifies the algorithm, since we 
do not have to distinguish between the initial and the 
further steps. The only difference with the genuine SF 
construction is that at time t one has mi + TOo/2 (instead 
of mt) edges present; hence, the asymptotic behavior of 
both models for t — > oo is the same. 

Three examples of the evolving networks of such a kind 
are given en Fig.l. Here is m = 3, L = 10 6 , N = 102 
and too = 6, (so that all three networks have exactly 300 
edges). Three different values of a were used: a = 0.0 
(scale-free model), a = 1.5 and a — 15.0. Note that in- 
creasing a leads to marked changes in the topology of the 
network. Fig. 1(a) corresponds to a genuine scale- free 
construction and exhibits a lot of long bonds connecting 
distance sites. On the other hand, only few such bonds 
are present in Fig. 1(c). 

In our further simulations we use a lattice of L = 2 ■ 10 7 
sites; the maximal number of the introduced nodes is 
N = 2 • 10 5 . All simulation results are based on the aver- 
age of 10 realizations of this structure. The error bars on 
Figs. 3-5 correspond just to this ensemble average. The 
simulations are done for several values of a and for two 
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values of m, the number of the outgoing bonds: m = 1 
and m = 3; mo = 2m. 
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FIG. 2. The degree distribution P(k) for different values 
of a and for m = 1 (a) and m — 3 (b). The values of a 
are a = (squares), a — 0.8 (crosses), a = 1.5 (triangles), 
a = 2 (filled circles), a = 5 (plusses) and a = 45 (diamonds). 
The dashed lines correspond to the theoretical curve for the 
scale-free model, (Ref. [1,12]) 



One of the prominent features of the scale free-model 
is that the distribution of the degrees of the nodes de- 
cays as a power law, i.e. like P(k) ~ fc -7 , with 7 = 3. 
This corresponds to the fact that the mean number of 
connections per site exists, but its dispersion diverges. 
Let us now discuss, how this distribution changes if the 
long-range connections are penalized. In Fig. 2 we plot 
the probability distribution of k for different values of a 
on double logarithmic scales. One readily infers that for 
all < a < 1 no important differences with the scale free 
model (a — 0) can be detected: in any case the asymp- 
totic behavior of P(k) is well-described by P{k) ~ k~ 3 . 



The distributions seem to be almost identical; however, 
small but statistically significant deviations can be de- 
tected for small k- values. At a ~ 1 the degree distri- 
bution shows a pronounced change in its behavior and 
ceases being a power law; now the behavior of the model 
with distance penalties is quite different. 

Let us concentrate on the case a > 1 and try to de- 
scribe the shape of the degree distribution under such 
conditions. The analysis of the simulations suggests that 
the corresponding mathematical expression could be a 
stretchcd-cxponcntial function of the form: 



P(k) = aexp(-6fc 7 ), 



(3) 



where the parameters a, b and 7 depend on a and m. To 
obtain the values of these parameter and to analyze the 
goodness of this fitting function we have fitted the data 
to Eq.(j^) using the nonlinear least-squares Levenberg- 
Marquardt algorithm p^| , taking into consideration the 
error bars as coming out of 10 realizations of each situ- 
ation. The data is replotted together with the outcomes 
of the fits in Fig. 4 on the scales in which the fitting 
function, Eq.(|^), is represented by a straight line. One 
namely takes fc 7 as the abscissa and lnP(fc) as the ordi- 
nata of the graph. Fig. 4 shows that such a fit (straight 
line) is surprisingly good! 

The values of the exponent 7 are shown as a function 
of a (a > 1) in Fig. 3, for the two different situations 
corresponding to m = 1 and m = 3. We see that 7 
monotonously grows with a, and that the dependences 
m = 1 and m = 3 differ, i.e. that the 7(a) dependence 
is nonuniversal. 




FIG. 3. The parameter 7 as a function of a. The upper 
dependence corresponds to m = 1, and the lower one m = 3. 
The lines are drawn as a guide for eyes. 



We note that in related models of growing networks 
another form of degree distribution appears: an expo- 
nentially demped power-law p3], 
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FIG. 4. Shown is lnP(fc) as a function of k 1 , where 7 is the output of the fit, Eq.(3). (See text for details). The parameters 
are: (a): m = 1, a = 1.5, 7 = 0.37. (b): m = 3, a — 1.5, 7 = 0.33. (c): m = 1, a = 3, 7 = 0.69. (d): m = 3, a = 3, 7 = 0.64. 
(e): m = 1, a = 10, 7 = 1.07. (f): m = 3, a = 10, 7 = 0.96. 



P(k) = afc 7 exp(-6fc). (4) 

We tested also this fit function and found out that it gives 
a good fit for larger a-values, but is definitely inferior to 
our fit, Eq.(g), for 1 < a < 3. 

A growing network with disadvantaged long bonds is a 



very interesting hierarchical construction. Thus, for large 
a, the strong correlation between the age of the connec- 
tion and its length exists. The old connections, made 
when the nodes were sparce, are typically long, while the 
younger connections get shorter and shorter, since more 
sites in the immediate vicinity of a newly introduced site 
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can be found. The simulations show that for a large, 
the nodes are almost surely connected to their nearest 
neighbors. On the other hand, the old, long-range con- 
nections are of great importance for the overall topology 
of the lattice, since they guarantee that for any a the 
network is a small- world one. 





ues of these labels divided by N — 1. This procedure is 
repeated for each node, and the overall mean value, the 
so-called path diameter of the network (I), is evaluated. 
The error bars of the figure correspond to the average of 
the mean diameters over 10 realizations of the network. 
Fig. 5 shows that the mean diameter of the network grows 
linearly in In TV, i.e. it shows the typical small- world be- 
havior. This behavior is preserved for all tested values 
of a; the largest value tested was a — 45, which, for 
m = 1, corresponds to a practically sure connection of 
a newly introduced node to its nearest neighbor. The 
high-a networks resemble closely the simple small- world 
constructions p5| . 

Let us summarize our findings. We considered a grow- 
ing network, whose growth algorithm is based, as in 
the SF construction, on a preferential attachment of the 
newly introduced nodes to the highly connected old ones. 
However, here the too long connections are disadvantaged 
by introducing penalties. Thus, the probability to con- 
nect two nodes separated by a distance d is proportional 
to d~ a , where a is a variable parameter. We found out 
that for a < 1 the degree distribution P{k) decays, as 
in the SF model, like P{k) ~ fc 3 , whereas for a > 1 a 
stretched exponential form P(k) — aexp(— bk 1 ) gives an 
extremely good description of this distribution. On the 
other hand, the small-world property is preserved at all 
checked values of a. 

The authors are grateful to Mr. F. Jasch for fruitful 
discussions and to Mr. M. Obermayr for the valuable 
technical assistence. Financial support by the Fonds der 
Chemischen Industrie is gratefully acknowledged. 
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FIG. 5. The diameter of a network as a function of the 
number of sites N. Panel (a) corresponds to a = 1.5 and panel 
(b) corresponds to a = 5. The upper lines in each panel are 
those for m = 1, the lower lines correspond to m = 3. Note 
the logarithmic scale. 



In Fig. 5 we plot the mean number of connections be- 
tween each two nodes of the network for two different 
values of a (a — 1.5 and a = 5) and for the two values 
to = 1 and m = 3 as a function of the network size N. 
The algorithm here is trivial: starting from a node (la- 
beled 0) we pass to all nodes connected to it (nodes of the 
first generation, labeled 1), then to nodes of the second 
generation (labeled 2), etc; untill all nodes are labeled. 
The mean distance between this node (labeled 0) and any 
other given node of the network is then the sum of all val- 
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